
subjectList = {...
    'K091412S'...
    'R091412V'...
    'M101012R'...
    };

numVox = 200;

extractedValues = zeros(1,6);


for subjIndx = 1:length(subjectList);
    
    [adaptationCurve] = stimHistory(subjectList{subjIndx}, numVox);
    
    display(adaptationCurve);
    
    figure(1);
    subplot(length(subjectList),1,subjIndx);
    ylabel('Average Beta (FOBS ROI)');
    for i=1:6
        extractedValues(i) = adaptationCurve{2,i+1};
    end
    plot(extractedValues);
    axis([1 6 -10 40]);
    hold on;
    
    figure(2);
    subplot(length(subjectList),1,subjIndx);
    ylabel('Average Beta (ADAPT ROI)');
    for i=1:6
        extractedValues(i) = adaptationCurve{3,i+1};
    end
    plot(extractedValues);
    axis([1 6 -10 40]);
    hold on;
    
    figure(3);
    subplot(length(subjectList),1,subjIndx);
    ylabel('Average Beta (Intersection ROI)');
    for i=1:6
        extractedValues(i) = adaptationCurve{4,i+1};
    end
    plot(extractedValues);
    axis([1 6 -10 40]);
    hold on;
    
end

figure(1);
subplot(3,1,1);
title('FOBS ROI');
subplot(3,1,3);
xlabel('Steps back');

figure(2);
subplot(3,1,1);
title('ADAPT ROI');
subplot(3,1,3);
xlabel('Steps back');

figure(3);
subplot(3,1,1);
title('Intersection ROI');
subplot(3,1,3);
xlabel('Steps back');
